function [vN1] = BoundaryPut(S, K, r, T, tau)

vN1 = exp(-r * (tau)) * max((K - S), 0);

